Entanglement Generated by the Dispersive Interaction: The Dressed Coherent State 
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In the dispersive regime of qubit-cavity coupling, classical cavity drive populates the cavity, but 
leaves the qubit state unaffected. However, the dispersive Hamiltonian is derived after both a frame 
transformation and an approximation. Therefore, to connect to external experimental devices, the 
inverse frame transformation from the dispersive frame back to the lab frame is necessary. In this 
work, we show that in the lab frame the system is best described by an entangled state known 
as the dressed coherent state, and thus even in the dispersive regime, entanglement is generated 
between the qubit and the cavity. Also, we show that further qubit evolution depends on both the 
amplitude and phase of the dressed coherent state, and use the dressed coherent state to calculate 
the measurement contrast of a recently developed dispersive readout protocol. 


The interaction between a two level system (TLS) 
and quantized electromagnetic radiation has been stud¬ 
ied extensively since the beginnings of quantum mechan¬ 
ics, with much effort devoted to the study of physi¬ 
cal systems described by the Jaynes-Gummings Hamil¬ 
tonian [1]. Over the last few decades the fields of cav¬ 
ity quantum electrodynamics (GQED) and more recently 
circuit quantum electrodynamics (cQED) have signifi¬ 
cantly developed, allowing for the exploration of the 
Jaynes-Gummings interaction in a wide range of parame¬ 
ter regimes and physical systems. In particular in cQED, 
both the strong coupling regime {g ^ k, 7, first achieved 
in Rydberg atoms [2]) and the strong dispersive regime 
(x At, 7) have been reached within the last decade [3]. 
In cQED, a superconducting qubit serves as the TLS, 
while the quantized electromagnetic fields are microwaves 
in either a strip-line resonator or 3D microwave cavity. 

Gontemporary experiments in cQED often work in the 
strong dispersive regime, where the qubit and microwave 
cavity are off resonance, and the Jaynes-Gummings inter¬ 
action reduces to an effective second order shift in system 
eigen-energies. In this regime, a wide range of quantum 
information protocols has been demonstrated [4], includ¬ 
ing quantum teleportation [5], entanglement generation 
by measurement and feedback [6, 7], non-classical mi¬ 
crowave state generation [8], and error correction by sta¬ 
bilization measurements [9]. 

When an empty electromagnetic cavity is driven by 
classical radiation, the state of the cavity is described 
quantum mechanically by the coherent state ja), where 
the complex amplitude a depends on the strength and 
length of the classical drive. In the dispersive regime 
of qubit-cavity coupling, when a classical cavity drive is 
applied the state of the joint system is typically described 
by the product state a \g) |ag) -I- b \e) |ae), with no qubit- 
cavity entanglement generated if the qubit is not initially 
in a superposition state. 

What is often overlooked is that the state a\g) \ag) + 
h\e)\ae) is an accurate description of the joint system 
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state under the dispersive approximation, which involves 
a frame transformation to the dispersive frame, and thus 
this state is not an accurate description of the qubit- 
cavity system in the lab frame of the experiment. In this 
manuscript, we will show that in the lab frame of the 
experiment a more accurate description of the joint state 
is the dressed coherent state \gje,a) [10]. Unlike the 
description in the dispersive frame, the dressed coherent 
state of the lab frame is entangled, even if the qubit is 
not initially in a superposition state. This has profound 
implications on the future evolution of the system, and 
we will show that future qubit evolution is dependent on 
both the amplitude and phase of the dressed coherent 
state. Similar effects have previously been studied for a 
driven qubit-cavity system where the qubit and cavity 
are resonant [11, 12[. 

This manuscript is organized as follows: in section I we 
describe the physical system of interest and the unitary 
frames we will be working in; in section II we define the 
dressed coherent state and give analytical and numerical 
evidence that this is the state created after the coupled 
system is driven classically through the cavity degrees 
of freedom; in section III we discuss applications of the 
dressed coherent state to quantum information protocols; 
finally, in section IV we make concluding remarks. 


I. THE PHYSICAL SYSTEM 

We consider a qubit and cavity coupled via the Jaynes- 
Gummings interaction, described by the lab frame Hamil¬ 
tonian 

H = Lojj^d - ' 2 ^^ ^ , ( 1 ) 

where d and are the usual bosonic annihilation and 
creation operators, is the Pauli matrix whose eigen¬ 
states are the qubit logical states, d^ are the qubit rais¬ 
ing and lowering operators, Wc/q are the cavity and qubit 
frequencies, g is the Jaynes-Gummings coupling strength, 
and we set h = 1 from here on. The eigenbasis for this 
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Hamiltonian is given by the dressed states [13] 


\g, n) = cos 0n \g, n) - sin 6 »„ |e, n - 1 ), ( 2 ) 

|e, n — 1) = cos9„ |e, n — 1) + sin0„ \g, n ), (3) 


to as the “dispersive-interaction” frame), after the rotat¬ 
ing wave approximation the full system Hamiltonian is 
then 


Hi(t) = -f (11) 


where the unbarred kets are the eigenstates of the un¬ 
coupled system, and the mixing angle is given by the 
relation 


On = - arctan (2X^/n ), (4) 

where A = 5 /A with A = Wq — Wc the cavity-qubit de¬ 
tuning. We work in the dispersive regime, defined by 
|A| <C 1. In this regime the mixing angle is well approx¬ 
imated as On ^ A-^n, provided n ^ ncriti where ricrit is 
the critical photon number at which the approximation 
breaks down due to the product \/nA approaching unity 
[14]. The dressed states then reduce to 

|g, n) = cos {X^/n) \g, n) — sin {X\/n) |e, n — 1), (5) 

|e, n — 1) = cos (Xy/n) |e, n — 1) -I- sin (Xy/n) \g, n) , (6) 

In this regime, it is then possible to transform to the 
dispersive frame by applying the unitary rotation Ud = 
exp IA (a^a — a^ai)}, and keeping terms up to first or¬ 
der in the dispersive shift x = 9^I The result of this 
procedure is the dispersive frame Hamiltonian 

ilD=Wca'a- - —(7) 


where 5 = Wc ~ Wd is the cavity-drive detuning, and 
the subscript “I” labels the interaction picture. In the 
dispersive-interaction frame, after a time T the state of 
the full system is described by 




f7i(T,0)|V>'(0)) 



IV''(0)), 


( 12 ) 


where 'T is the usual time ordering operator. 

To calculate Ui{T,0) exactly we will use the Magnus 
expansion [15], given by 


Ui{T, 0) = exp 


J2^k{T,0) 


\k—l 


(13) 


where Xlk{T,0) is the /c’th order Magnus generator. For 
our system these generators are zero for k > 2, and so 
tJi{T,Q) can be calculated exactly. The first order Mag¬ 
nus generator is given by an expression proportional to 
the average Hamiltonian 

ni(T, 0) = -i [ dt Hi{f) = \g) ( 5 I (ag(T)a'l' - ag(T)a) 
Jo 

+ \e) {e\{ae{T)d^ - al{T)d), (14) 


which we emphasize is not in the lab frame, but in the 
dispersive frame defined by Ud. To highlight the in¬ 
equivalence of the two frames, we note the identities 
(which will be useful later) 

Ul\g,n) = \g,n), ( 8 ) 

Ul\e,n) = \e,n), (9) 

where the dressed states are given by equations (5) and 

( 6 ). 


II. THE DRESSED COHERENT STATE 
A. Analytic Derivation 


where 


cxg{T) 


ae{T) 


-e* - 1 ) 

S-X 

-e* - 1 ) 

S + x 


(15) 


As [Hi{ti), Hi{t 2 )] oc /(cTz) 0 1 which commutes with 
the Magnus expansion truncates at second order, 
and the second order generator results in a qubit-state 
dependent relative Stark phase. Also, as 142 (T,0) com¬ 
mutes with Oi(T, 0 ) the full evolution operator is 


i7i(r, 0) = exp {F!i(T, 0)} exp 0)} 

= (iff) {g\ b{ag{T)) + \e) (e| b{a,{T))) (16) 


We are interested in the effect on the full system of a 
classical drive applied to the cavity, as described in the 
lab frame and to lowest order in the dispersive frame by 
the Hamiltonian 

= 2 cos(u;dt) (ea -f e*o^) , ( 10 ) 

where Wd and e are the frequency and complex amplitude 
of the drive respectively. In the interaction frame of the 
system Hamiltonian Hd of equation (7) (which we refer 


where £>{(}) = exp {/3a^ — /3*a} is the usual displacement 
operator, and is the phase from the fl 2 iT, 0 ) 

term, shown explicitly in equation (A2). 

In the lab frame of equation (1) we consider an ini¬ 
tial state given by [ 5 , 0 ), which is the ground state of 
the Jaynes-Cummings Hamiltonian. This state is a dark 
state and therefore invariant under both the transfor¬ 
mation into the dispersive frame UD(t) and the trans¬ 
formation into the interaction frame exp liHDt\, and 
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so in the dispersive-interaction frame the initial state is 
|'0'(O)) = |5,0) (we work in the lab frame basis through¬ 
out). 

We can then calculate the final state in the dispersive- 
interaction frame 

IV'D(r)) = (|5) (5l D{a^{T)) + |e) (e| b{a,{T))) \g, 0) 

= \g,oc^{T)), (17) 


where as our initial state is not a qubit superposition 
state the effect of the second order Magnus term is a 
global phase that can be ignored. Next we transform 
this trivially out of the interaction frame back to the 
dispersive frame 


|^D(r)) = exp{-*ffDT}|V-()(T)) 




g,ag(r)e-*(--^)^\ = |5,5g(T)). (18) 


To return to the lab frame we must apply the inverse 
dispersive transformation and using equation (8) we 
see that the final state in the lab frame is 


IV-(r)) = Ul |5,5g(T)) = ^ \9,n) 

„ vn! 

(19) 





Figure 1: Schematic diagrams of the dressed coherent states 
\g,a) and |e, a) for \a\^ = 4. The red and blue curves are 
the Poisson dis tribution of superposition coefficients, = 
K 5 /e,n| g/e,a)\^. 


B. Effective Qubit Drive 


The qubit-cavity state of equation (19) is the dressed 
coherent state \g,ag{T)), where in general a dressed co¬ 
herent state has the form 


To understand the dressed coherent state for A 1 it 
is useful to expand the state of equation (20) in powers 
of A to obtain 


_ I ,2 _ fyll' _ 

\g/e, a) = e ^ ^ ^ l^/e, n) 

„ vn! 


( 20 ) 


1-oV -IW P 


(l - \9,n) 


where \g/e,n) is given exactly be equation (2)/(3), and 
to first order in A by equation (5)/(6). To the best of 
our knowledge this state was first described in Ref. [10]. 
Schematic diagrams of the dressed coherent states \g, a) 
and |e, a) for |ap = 4 are shown in FIG. 1. The curves 
in FIG. 1 represent the Poissonian weights of the sum 
in equation (20), and highlight the fact that the dressed 
coherent state has the same distribution of superposition 
coefficients as the coherent state, only for the dressed 
states instead of the bare states. 

As we have just shown, if the system starts in the state 
\g, 0) then the final state of the system after a classical 
cavity driv e of length T will be the dressed coherent state 
|g, Qfg(T)). Similarly, if the system starts in its first ex¬ 
cited state, given by |e, 0) in the lab frame, then the final 
state after classical drive will be the state |e, de(T)) (see 
appendix A for further details). We will now discuss an 
intuitive physical understanding of the dressed coherent 
state, and compare the analytic results to full numerical 
simulations. 


— X\Pn |e, n — 1) 


+ O (A^) 


= (Iff)-A/3 |e)) 1/3) 

( 21 ) 


and we see that to lowest nontrivial order in A the qubit is 
effectively in the superposition state ([(;) — A/3 |e)) /\fM, 
with the normalization factor N = 1 -I- A^|/3p. Thus, to 
lowest order in A, when the qubit-cavity system is driven 
by a classical cavity drive (creating a dressed coherent 
state) the qubit is effectively weakly driven on resonance 
via the Hamiltonian (in the interaction picture) 

^Eff = goo-” + gSd’+, (22) 


for a time T as before, where the effective qubit drive 
strength is given by go = {i/3*X)/T for the dressed coher¬ 
ent state \g, /3). 
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It is important to note that under the dispersive ap¬ 
proximation the cavity drive of equation (10) is only the 
leading order term, and that the next order correction in 
the dispersive frame, of order A, is a qubit drive. How¬ 
ever, this drive is at a frequency Wd, which will not be res¬ 
onant with the qubit transition frequency, and as a result 
this qubit drive has no net effect on relevant timescales. 
For this reason we have not included this off-resonant 
qubit drive in our analtyical calculations. 

The effective resonant qubit drive described here by 
equation (22) is not due to this first order term in the 
dispersive approximation of the cavity drive Hamiltonian, 
as it occurs even when only the zeroth order term of the 
dispersive frame cavity drive is considered, as in equation 
(10). It is uniquely an effect of considering the lab frame 
state for qubit-cavity interaction with a driven cavity, 
and results from interactions between the cavity and the 
qubit. 

C. Numerical Simulations 

In section HA we have shown that under the disper¬ 
sive approximation done completely, the final state in the 
lab frame of a qubit-cavity system after a classical cavity 
drive is a dressed coherent state. While for A ^ 1 and 
|ap ^ Ucrit it is sufficient to keep terms up to order A^ 
and obtain the Hamiltonian of equation (7), it is worth¬ 
while to examine what effect the neglected higher order 
terms have on the final state. To do so, we numerically 
investigate time evolution induced by the Hamiltonian 

Hrit) = WcO^Q,- + 9 {d'~a^ + a'^a) 

+ (ee“‘'*a -f 0(t - T), (23) 


figures. Fidelities for both the ground or excited dressed 
coherent state are plotted. Due to decreasing validity of 
the dispersive approximation, both FIGs. 2(a) and 2(c) 
show a decreasing overlap as either \a\’^ or A increase. As 
is to be expected, increasing A is more detrimental to the 
agreement between the numerical state and the dressed 
coherent state, as a larger A requires fewer photons in the 
cavity for terms beyond order A^ in the full Hamiltonian 
to become relevant. For increasing jap, FIG. 2(a) shows 
that even for a high average photon number of jap = 9 
the numerical state still has an overlap greater than 90% 
with a dressed coherent state. 

Interestingly, FIG. 2(d) shows that for increasing drive 
strength e, the overlap between the dressed coherent state 
and the numerical state actually increases. This can 
likely be understood by competition between the cav¬ 
ity drive and higher order effects beyond the dispersive 
approximation, which will both be off-diagonal in the ba¬ 
sis of equation (7). The stronger the drive, the more it 
dominates this competition, which therefore diminishes 
the effect of the higher order correction terms, leading to 
a state closer to the dressed coherent state. 

Finally, FIG. 2(b) shows the difference in fidelity be¬ 
tween the dressed coherent state and the undressed prod¬ 
uct state, i.e. e. A) — e, A), as a func¬ 

tion of |ap. As can be seen, the difference is always 
positive, and the dressed coherent state is always a bet¬ 
ter description of the numerical state than the undressed 
product state. 

III. IMPLICATIONS AND APPLICATIONS 
A. Qubit Operations 


for the lab frame initial states |( 7 , 0 ) and |e, 0), with 
Wd = Wc ~ X Wd = Wc -f X respectively. We simu¬ 
late over a range of A, e, and target dg/e(7’) to see how 
these parameters affect the accuracy of the dressed co¬ 
herent state. 

We compare the numerically created states with the 
dressed coherent states |g, dg(T)) and \e,ae{T)), as well 
as the undressed product states |( 7 ,dg(T)) and |e, Q;e(7")) 
(all in the lab frame), by calculating the overlaps 


■^D^''(|aP.e,A) 

.Fs/^(|a|2,e,A) 


(V’(T)I |5/e,dg/e(T)) 
(V’(T)I |5/e,dg/e(T)) 


(24) 

(25) 


where \'>p(T)) is the state created by numerical simula¬ 
tion, and we have set that for either initial state the tar¬ 
get coherent state amplitude is the same, i.e. |Q;g(r)| = 
|de(T)| = |q;(T)|. The phase of a{T) has no impact on 
the fidelity and so the fidelity depends only on |ap. 

Figures 2(a), 2(c), and 2(d) show 1 — {\a\'^, e, X) for 

varying |ap. A, and e respectively, with the non-varying 
parameters held at the constant values indicted on the 


In this section we examine a simple case of qubit op¬ 
erations performed on a dressed coherent state to em¬ 
phasize the dependence of these operations on both the 
amplitude and phase of the cavity state. Such interac¬ 
tions could occur in a set-up that involved parallel oper¬ 
ations on a qubit and a cavity used as a bus [16] or as a 
quantum memory [17]. We consider a primitive for qubit 
operations that, while itself not necessarily useful or in¬ 
teresting, elucidates the relevant physics involved, such 
that useful and interesting applications may be developed 
from it in future work. Starting in the dressed coherent 
state \g,P) (created as described in section HA), we at¬ 
tempt to rotate the reduced state of the qubit to as close 
to the excited state as possible. We will show that the 
excitation probability will depend on both the amplitude 
and phase of /3. 

We begin by expanding the inital dressed coherent 
state |g,/3) to lowest order in A, as in equation (21). In 
this case, the reduced state of the qubit is a weakly ro¬ 
tated qubit state, with the phase of this rotation propor¬ 
tional to the phase of /3. For a weakly rotated qubit state, 
the probability amplitude upon further rotating the qubit 
depends on the phase of the applied drive relative to the 
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Figure 2: Comparison of the numerical state created by simulation of equation (23) with the dressed coherent state as a function 
of (a) photon number |ap, with e and A constant, (c) A, with e and |ap constant, and (d) drive strength e, with \a\^ and A 
constant, (b) The difference in overlap with the numerical state between the dressed coherent state and the undressed product 
state as a function of photon number |ap, with e and A constant. Values of \a\‘^, A, and e chosen to be commensurate with 
current experiments. Additional phase matching was required to obtain high fidelity (see appendix B). 


phase of the initial qubit rotation, and therefore, the suc¬ 
cess of rotating to the excited state from the state \g, f3) 
will depend on the phase of /3. For purely real (imagi¬ 
nary) /3, the reduced qubit state lies in the x-z (y-z) plane 
on the Bloch sphere, and rotation about the x (y) axis 
cannot transform the qubit state to the state |e), while 
rotation about the y (x) axis can. For general /? the su¬ 
perposition is along the plane defined by the z axis and 
the line {{x,y) = {cos{(j)p)s,sin{(j)p)s), Vs G K} in the x- 
y plane, which depends on the phase of /3, given by (fip. In 
this case perfect state transfer to |e) is not possible by ro¬ 
tation along any axis other than that defined by (pp+i:/2, 
given by the line {{x,y) = (sin((/)^)s, cos(0^)s). Vs G K}. 

More rigorously, we consider the qubit-drive Hamilto¬ 
nian to first order in the dispersive frame given by 

ijQ = ,7e-“‘<T+ -b . (26) 

In the interaction picture with respect to the dispersive 
Hamiltonian of equation (7), this Hamiltonian becomes 

H'^ = -b , (27) 

where fi = a) a and v = ojq + x ~ To calculate the 
evolution operator, we follow the same procedure as in 
section HA using the Magnus expansion. Unfortunately, 


the Magnus generators do not vanish for any finite 
order n, due to the nonlinearity of the qubit. To account 
for this we choose our evolution time to be short enough 
(roughly one period for the drive frequency u) that Mag¬ 
nus terms beyond first order have minimal contribution 
to the evolution, as the magnitude of their effects is small 
and only becomes relevant after accumulating for a long 
period of time. With this in mind, the analytical results 
presented below are meant to be instructive and to high¬ 
light the important physical effects, rather than to have 
high accuracy. 

The result of this calculation (see appendix C for fur¬ 
ther details) is a very complicated expression for the ex¬ 
cited state probability as a function of time, given by 
equation (C5). In order to gain some intuitive under¬ 
standing of the rigorous result, we set uj = cjq and go to 
the limit where x, A —> 0, in which case, the probability 
of finding the qubit in the excited state at time r is given 

by 

Pe('r) = (l — sin^ ( tjt ) + A^|/3|^ cos { 2 t ] t ) 

-b A sin(2ryr) |^lm(,0e“*‘^) cos(Ar) 

-b Re(/3e“*‘^) sin(Ar)^, (28) 
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Figure 3: Qubit excited state probability as a function of the 
time of the applied qubit drive, starting in a dressed coherent 
state with either purely real or purely imaginary /? for |/3p = 
4. The drive strength is strong, It;! = O.OSoiq, so that phase 
effects are clearly visible. The drive frequency is set to the 
average shifted qubit frequency, uj = tJq + 2(|/5|^ + l), rj is 
purely real (i.e. \r]\ — rj), and A = 0.1 as elsewhere. 


where = r]/\r]\, and A = Wq — Wc as before. From the 
last term in equation (28) we see that the excited state 
probability depends not only on the photon number in 
the cavity, but also the interplay between the coherent 
state phase and the phase of the drive, through the un¬ 
equal dependence on Re(,de“*‘^) and Im(/3e“*‘^). This 
effect can be understood as interference between the ef¬ 
fective qubit drive of equation (22) caused by the dressed 
coherent state and the applied qubit drive of equation 
(26), in a manner similar to coherent destruction of tun¬ 
neling [18]. 

Unfortunately, while easy to understand, equation (28) 
is not very accurate in the relevant parameter regimes 
(due to the approximations made), and the error in trun¬ 
cating the Magnus expansion at first order grows for 
longer times. However, the analytical calculations lead¬ 
ing to equation (28) were done to distill the relevant 
physical effects and present them in an understandable 
manner, not to obtain highly accurate results. To accu¬ 
rately test the phase dependence effects, we numerically 
simulate a qubit-cavity system with classical qubit drive 
(without making the dispersive approximation). As an 
initial state, we start with the state created by simulation 
of equation (23) of section IIC, which has high overlap 
with a dressed coherent state. Figure 3 shows the ex¬ 
cited state probability as a function of time, starting in 
a dressed coherent state with either purely real or imag¬ 
inary /3, and with a purely real qubit drive rj. 

In addition to photon number effects, where the qubit 
resonance frequency is modified by the photon number 
in the cavity [19], FIG. 3 shows that the excited state 
probability’s time evolution depends also on the phase 
of the dressed coherent state amplitude /3. This agrees 
with the intuitive conclusions drawn previously, using the 
approximation of equation (21), where one considers the 
dressed coherent state to be a coherent state in the cavity, 
and a weakly rotated qubit. 


B. Dispersive Multi-Qubit Readout with a 
Threshold Detector - Limited Contrast 

In recent proposals for dispersive single qubit readout 
[20], and qubit parity readout [21], protocols were devel¬ 
oped that conditionally populate the cavity dependent 
on the state (parity) of the qubit(s), after which, by us¬ 
ing a threshold photon counter to distinguish between 
the qubit dependent cavity states, the state (parity) of 
the qubit(s) can be non-destructively measured. A sub¬ 
sequent coherent cavity drive of opposite phase removes 
the cavity occupation. Such conditional cavity occupa¬ 
tion can be achieved by setting Wd = Wc — x in equation 
(15), and the correct choice of T such that aAT) ^ 0 
while ae{T) = 0. 

For single-qubit readout, starting from an initial state 
in the lab frame of either jg, 0) or je, 0) and using equa¬ 
tions (19) and (A6), after classical cavity drive the qubit 
state dependent qubit-cavity states are 

IV-(T)) = \g,ag{T)), (29) 

|vl/(r)) = cos (A) ^ sin (A) \g) |C(T)), 

(30) 

for the suitable choice of T that ensures ae{T) = 0. 
Misidentification of the qubit excited state as the ground 
state occurs when there is spurious photon population in 
the cavity when the qubit is in its excited state. As can 
be calculated using equation (30) this spurious photon 
population is given by 

N = (4'(T)| a^a |«'(r)) « sin^ (A) cos^ (A) (1| a^a |1) 

+ sin2 (A) {g\{aT)\Uua^aUl\g)\aT)) 

« sin^ (A) (cos^ (A) -f (C(T)| al'a \^(T))) 

= sin^ (A) (cos^ (A)-I-1-f jagp) (31) 

where in the first approximation we have ignored the 
cross terms as the photon occupation of |C(r)) is much 
greater than 1, and in the second approximation we have 
assumed the dispersive transformation only slightly mod¬ 
ifies the average photon number of |^(T)). The photon 
occupation of |C(r)) is high since it is the single photon 
Fock state displaced by D{ag{T)), and ag{T) > 0 for this 
protocol. Extending this effect to multiple qubits in their 
excited states explains the even parity misidentification 
error for four qubit parity measurement seen in Ref. [21], 
as it is the four-qubit generalization of the state |‘C(7^)) 
that leads to spurious detections when the qubits are in 
an even parity state. The exact value of this error is de¬ 
pendent on the nature of the threshold photon counter 
used [22, 23]. 

IV. CONCLUSION 

In the work presented here, we have undertaken a crit¬ 
ical and careful examination of the dispersive limit, dis¬ 
persive approximation, and the dispersive Hamiltonian, 
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in order to understand the distinction and commonalities 
between these commonly used terms. We have found that 
one obtains the dispersive Hamiltonian after making the 
dispersive approximation, which requires both being in 
the dispersive limit of the Jaynes-Cummings interaction, 
and a frame transformation to the dispersive frame. The 
dispersive Hamiltonian therefore describes evolution of 
the system in the dispersive frame, and in order to ob¬ 
tain results valid in the lab frame, one must apply the 
inverse frame transformation. 

This has a profound effect on the description of the 
system state, and in particular, we have found that after 
a classical drive is applied to the cavity, the state of the 
system in the lab frame is accurately described by the 
dressed coherent state \g/e,a), not, as often used, the 
product state \g/e) ja). The entanglement present in the 
dressed coherent state will affect all future operations 
on the qubit and the cavity. We have shown how this 
is relevant to rotations of the qubit state, in particular 
that the probability of rotating the qubit from its ground 
to excited state (or vice versa) depends on the phase 
difference between a and the applied qubit drive. We 
have also explained the limit in measurement contrast for 
readout of the qubit state via a cavity and a threshold 
photon counter, reported in Ref. [21]. Future work will 
continue to explore the effects of the dressed coherent 
state on contemporary quantum information protocols. 
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Appendix A: Dressed Coherent State for An Excited 

Qubit 

In this appendix we consider starting in the initial state 
in the lab frame of the bare excited qubit state |e, 0 ), 
which can be prepared by initialization to the state \g, 0 ) 
followed by a fast, nonadiabatic pulse on the qubit. In 
the dispersive frame the initial state will be |4>d(0)) = 
[/ole, 0) = cos (A) I e, 0)— sin (A) |( 7 ,1). Transforming into 
the interaction frame, the state is unchanged since t = 0. 
To calculate the final state, we will need the second order 
Magnus generator, given by 

^^2(T,0) = -i^ dh 

= (Al) 

From this we can calculate the qubit-state dependent 
phase due to the second order Magnus generator, given 


dU 




by 


^^F(a,T) ^ ( sin(xT(j^) - xTctz)^ , (A2) 

and make the identification 

F{azT) = ^(^sin(xTcr2) - (A3) 

Now using equation (16), at t = T the system state in 
the dispersive-interaction frame is given by 

|vI/'D(T))=cos(A)|e,ae(T)) 

_ edF(T)-Fi-T)] 1^^ D{a^iT)) |1), (A4) 

where we have factored out a global phase. Transforming 
back into the dispersive frame we have: 

I^D(r)) =cos(A)|e,ae(T)) 

- (T) e-*(‘"“-^"^)“'“^.D(ag(r)) \g, 1) 

= cos (A) |e, <5e(r)) - sin (A) \g) \^{T)) , (A5) 

where de(P) = ae(T)e-*(‘^-+x)^, G(T) = (wq -k x)T + 

2F(T), and |e(T)) = |1) is 

the displaced single photon Fock state. Finally, trans¬ 
forming back into the lab frame gives 

|4<(r)) = ul (cos (A) |e, de{T)) - sin (A) \g) |5(t))) 

= cos (A) |e, <5e(T)) - sin (A) \g) , (A6) 

where as before we have made the identification that 
Uq |e, de(T)) = |e, de(T)) to first order in A. If A ^ 1 
such that sin(A) « 0, then as before |4'(T)) is a dressed 
coherent state. Additionally, if as an initial state we use 
|4'(0)) = |e, 0) instead of the bare excited state |e, 0), 
then as Uj:,\e,0) = |e, 0) the term proportional to sin(A) 
in equation (A6) disappears, and |4>(r)) contains only 
the dressed coherent state |e, de(T))- The dressed initial 
state is closer to experimental reality as it is an eigenstate 
of the Hamiltonian, however, the calculation for the bare 
initial state |'I'(0)) = |e, 0) was shown for completeness, 
as some experimental protocols can prepare this state, 
and as from the solution for the bare initial state the 
solution for the dressed initial state is trivial to obtain. 



Appendix B: Corrections to the Phase of the 

Dressed Coherent State Due to Nonlinear Terms 

In the expected parameter regime, FIG. 2 demon¬ 
strates excellent agreement between the state created by 
a numerical simulation of the Jaynes-Cummings Hamil¬ 
tonian and the dressed coherent state. However, to 
obtain this high fidelity, it was necessary to include 
effects beyond the dispersive Hamiltonian to correctly 
match the phase of the dressed coherent state amplitude 
Q!g/e(T) with the numerical state. In particular, follow¬ 
ing Ref. [24], we include the nonlinear term proportional 
to AA^ in the classical equations of motion for the am¬ 
plitudes ag/e(T) and obtain the approximate solutions 

ag(r) = ag(T) exp - x -f , (Bl) 

de(T) = Q!e(r)exp + X- C -f T^, 

(B2) 

where ^ = AA^, {n{T)) = (a’I’a) (T) is the average photon 
number in the cavity at time T, and ag/e{T) are defined 
as before in equation (15). Using these modified coherent 
state amplitudes in equations (19) and (A6) gives the ex¬ 
cellent overlap with the numerical states seen in FIG. 2. 
We emphasize that the nonlinear corrections have only 
been used to correct the classical solution for ag/e{T), 
and that to lowest order the nonlinearity modifies only 
the phase and not the amplitude of the dressed coherent 
states. Furthermore, the squeezing Hamiltonian of the 
nonlinearity is not considered in our analytical calcula¬ 
tions (as this would modify the state so that it was no 
longer a dressed coherent state), and while its effect is 
small, we suspect it to be the leading cause of the less 
than unit fidelity seen in FIG. 2. 


Appendix C: Driven Qubit Excited State Probability 

To first order in the Magnus expansion, the qubit-drive 
evolution operator is 

Uq{t,0) = exp{-iOi(T,0)} 

= cos (^\r]b{n, t)\{v + 2xh)~^^ I 

-f sin {\'qh{n, t)\ + 2xh)"^^ T))*cr" 

- T]b{n,T)a+'^\r]b{h,T)\~'^, (Gl) 


where we have defined the operator function 6(h, r) = 
1 — Starting from the initial state |V’d(0)) = 

\g,l3) in the dispersive-interaction frame, the state after 
a time t of the applied qubit-drive is 


IV’D('r)) = Uq{t,0) |^/’d(0)) 


-\0r 
= e 2 


E 




— I sin 


f \'nKk,T)\ 
E -h 2kx 






(G2) 


where = ri/\r]\, and we have used the fact that 
&(fc, T)/|6(fc, r)| = , Transforming out of the 

interaction frame, we arrive at the state in the dispersive 
frame (after factoring out a global phase) 


IV'D(r)) = U/(r,0) IV’D('r)) 


= e 


cos 




— I sm 


\v + 2kx J 


(G3) 


where (5 = absorbs the rotation of the cavity 

state in phase space, and S = Wq-|-x + w-|- 2kx = v + 
2kx + 2(jj. Finally, transforming back into the lab frame, 
we end up with the state 


iV'D(r)) = uI){t,0) IV’d(t)) = 




T — 


COS 




— cos 


\i, + 2kxj V y^ 


— I sm 


— i sm 






1 ) 


(G4) 


Now we can calculate the probability at a given time r 
that the qubit is in the excited state, given by Pe(T) = 
(V'd(t)I |e) (e| Gl |V'd(t)). Using equation (G4) we calcu¬ 
late this to be 
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Equation (C5) is quite cumbersome, and to gain some x, A —)■ 0 limit, which results in equation (28). 

intuitive understanding, we set a; = Wq and examine the 


[1] E. Jaynes and F. Cummings, IEEE Proc. 51, 89 (1963). 

[2] D. Meschede, H. Walther, and G. Muller, Phys. Rev. 
Lett. 54, 551 (1985). 

[3] A. Wallraff, D. Schuster, A. Blais, L. Frunzio, R. Huang, 
J. Majer, S. Kumar, S. Girvin, and R. Schoelkopf, Nature 
431, 162 (2004). 

[4] M. H. Devoret and R. J. Schoelkopf, Science 339, 1169 
(2013). 

[5] L. Steffen, Y. Salathe, M. Oppliger, P. Kurpiers, M. Baur, 
C. Lang, G. Eichler, G. Puebla-Hellmann, A. Fedorov, 
and A. Wallraff, Nature 500, 319 (2013). 

[6] D. Riste, M. Dukalski, C. A. Watson, G. de Lange, 
M. J. Tiggelman, Y. M. Blanter, K. W. Lehnert, R. N. 
Schouten, and L. DiCarlo, Nature 502, 350 (2013). 

[7] J. M. Chow, J. M. Cambetta, E. Magesan, D. W. Abra¬ 
ham, A. W. Cross, B. R. Johnson, N. A. Masluk, C. A. 
Ryan, J. A. Smolin, S. J. Srinivasan, et ah, Nat Commun 
5, 4015 (2014). 

[8] G. Kirchmair, B. Vlastakis, Z. Leghtas, S. E. Nigg, 
H. Paik, E. Ginossar, M. Mirrahimi, L. Frunzio, S. M. 
Girvin, and R. J. Schoelkopf, Nature 495, 205 (2013). 

[9] J. Kelly, R. Barends, A. G. Fowler, A. Megrant, E. Jef¬ 
frey, T. C. White, D. Sank, J. Y. Mutus, B. Campbell, 
Y. Chen, et ah. Nature 519, 66 (2015). 

[10] E. A. Sete, A. Galiautdinov, E. Mlinar, J. M. Martinis, 
and A. N. Korotkov, Phys. Rev. Lett. 110, 210501 (2013). 

[11] V. Peano and M. Thorwart, EPL (Europhysics Letters) 
89, 17008 (2010). 

[12] V. Peano and M. Thorwart, Phys. Rev. B 82, 155129 

( 2010 ). 


[13] S. Haroche and J.-M. Raimond, Exploring the Quantum: 
Atoms, Cavities, and Photons (Oxford University Press, 
Oxford, 2006). 

[14] A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin, and 
R. J. Schoelkopf, Phys. Rev. A 69, 062320 (2004). 

[15] W. Magnus, Communications on Pure and Applied 
Mathematics 7, 649 (1954). 

[16] A. Galiautdinov, A. N. Korotkov, and J. M. Martinis, 
Phys. Rev. A 85, 042321 (2012). 

[17] Z. Leghtas, G. Kirchmair, B. Vlastakis, R. J. Schoelkopf, 
M. H. Devoret, and M. Mirrahimi, Phys. Rev. Lett. Ill, 
120501 (2013). 

[18] F. Grossmann, T. Dittrich, P. Jung, and P. Hanggi, Phys. 
Rev. Lett. 67, 516 (1991). 

[19] D. 1. Schuster, A. A. Houck, J. A. Schreier, A. Wallraff, 
J. M. Gambetta, A. Blais, L. Frunzio, J. Majer, B. John¬ 
son, M. H. Devoret, et ah, Nature 445, 515 (2007). 

[20] L. C. G. Govia, E. J. Pritchett, C. Xu, B. L. T. Plourde, 
M. G. Vavilov, F. K. Wilhelm, and R. McDermott, Phys. 
Rev. A 90, 062307 (2014). 

[21] L. C. G. Govia, E. J. Pritchett, B. L. T. Plourde, M. G. 
Vavilov, R. McDermott, and F. K. Wilhelm (2015), 
arXiv: 1502.03340. 

[22] L. C. G. Govia, E. J. Pritchett, S. T. Merkel, D. Pineau, 
and F. K. Wilhelm, Phys. Rev. A 86, 032311 (2012). 

[23] L. C. G. Govia, E. J. Pritchett, and F. K. Wilhelm, New 
Journal of Physics 16, 045011 (2014). 

[24] M. Boissonneault, J. M. Gambetta, and A. Blais, Phys. 
Rev. A 79, 013819 (2009). 



